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Abstract 

The anomalous dimensions of trilinear-quark operators are calculated at 
leading twist t = 3 by diagonalizing the one-gluon exchange kernel of the 
renormalization-group type evolution equation for the nucleon distribution 
amplitude. This is done within a symmetrized basis of Appell polynomials 
of maximum degree M for M 3> 1 (up to order 400) by combining analyti- 
cal and numerical algorithms. The calculated anomalous dimensions form a 
degenerate system, whose upper envelope shows asymptotically logarithmic 
behavior. 
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The anomalous dimensions of field operators are crucial ingredients in determining the 
scale dependence of hadronic wave functions in QCD. Though the calculation of the bound- 
state wave function at the initial (low-momentum) scale cannot be determined within per- 
turbative QCD and requires nonperturbative methods (for a review, see, [0), its evolution 
is governed by a renormalizat ion-group type evolution equation which in leading order 
reduces to the one-gluon-exchange kernel. 

Bilinear quark operators of lowest twist (t = 2) are representations of the coUinear con- 
formal group and therefore their anomalous dimensions (to first order in ag) are completely 
determined by collinear conformal invariance (see, e.g., Therefore, the eigen- 

value problem of twist-2 quark operators with an arbitrary number of total derivatives is 
exactly solvable. In contrast, the anomalous-dimensions matrix of trilinear quark opera- 
tors of leading twist-3 cannot be diagonalized completely by collinear conformal invariance 
|7|j8|J^,p!0|, p!l| , p!^ . As a result, neither the eigenfunctions nor the eigenvalues (alias the anoma- 
lous dimensions) of the nucleon evolution equation are known explicitly. Hence one has to 
diagonalize the anomalous dimensions matrix order by order within an appropriate operator 
basis. 

While operators involving total derivatives do not contribute to forward scattering am- 
plitudes, like deep-inelastic scattering, exclusive processes receive contributions from non- 
forward matrix elements, and consequently the mixing with trilinear quark operators of 
twist-3 containing total derivatives cannot be neglected. Given that the baryon (hadron) 
distribution amplitudes are determined at low normalization scales, of the order of 1 GeV 
or less, |ll| Jl3|JI^Jl5[| , the knowledge of anomalous dimensions is indispensable in confronting 
calculated observables, like the proton form factor, with experimental data at larger values 
of the momentum transfer. 

In this letter we present a systematic treatment of the nucleon evolution equation, fo- 
cusing on the large-order, i.e., quasi-asymptotic, behavior of the spectrum of the anomalous 
dimensions of twist-3 multiplicatively renormalizable baryonic /1/2-operators with the num- 
ber M of total derivatives acting on these fields tending to infinity. More precisely, we 
consider operators of the following form 

^Kn.n3)(Q) ^ . ^)ni ^(q) ^ (^^ . j^p ^(q) . j^^n, ^^^^^q^ 

with rii + n2 + = M, and where z is an auxiliary lightlike vector (i.e., z"^ = 0) to project 
out the leading twist contribution Note that since there are M + 1 independent operators 
at each order, one has to diagonalize a (M -|- 1) x (M + 1) matrix in order to determine 
the corresponding eigenvectors and eigenvalues. This is done here within a basis of trilinear 
quark operators in terms of appropriately symmetrized Appell polynomials (see below) by 
analytically diagonalizing the nucleon evolution equation up to order 7, and proceeding 
beyond that order up to a maximum order of 400 numerically. 

We have already elaborated on these issues in a series of papers [16,17, IS, 1^,^,1^, and 



the present work complements and generalizes our previous results, establishing in particular 
that the upper envelope of the spectrum for very large order M follows a logarithmic law 
given by the fit 

7„(M) = -0.638 + 0.8882 ln(M + 1.91) (2) 
for any member of the spectrum labeled by the index n. 
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Before proceeding with the presentation of our analysis, let us first mention a very re- 
cent alternative approach by Braun and collaborators |21], based on the observation that 
the eigenvalues of a helicity A = 3/2 baryon, (say, the A++) can be determined exactly 
because the corresponding Hamiltonian H3/2 is integrable and can be formally treated as 
a one-dimensional XXX Heisenberg spin magnet of noncompact spin s = — 1. Though the 
Hamiltonian 7ii/2 (relevant for the nucleon) is not integrable and consequently the corre- 
sponding eigenvalue problem is not exactly sovable, these authors argue that the helicity-flip 
terms can be treated as a small perturbation, so that the spectrum of anomalous dimensions 
can be determined by calculating the mixing matrix elements of these terms in the basis 
of the exact eigenf unctions of 7^3/2- This is an attractive mathematical framework, and we 
shall compare our results with theirs at the end of the paper. 

Let us now enter our approach in more detail. 

The evolution equation is given by 



X1X2X3 



where 



2Po 

7- / T2-"s(^± 



Po Jo 



[dy]V{x,,y,)^{y,,Q') , (3) 



In 



as(/i^ 



In 



InQVA^CD 



(4) 



is the evolution "time" parameter, /3o = 11 — |nf = 9 for three flavors, $ = X1X2X3 $ denotes 
the full nucleon distribution amplitude, {$} its eigenf unctions, and Cp = (A^"^ — 1) /2Nc, 
with the integration measure defined by [dy] = dyidy2dy36 (1 — |/i — ^2 — Us)- The kernel 
V{xi,yi) = V{yi,Xi) is the sum over one-gluon interactions between quark pairs {i,j} at 
order and is not a function but a distribution from which infrared divergences at Xi = yi 
have been removed. 

We propose to solve this evolution equation by employing factorization of the dependence 
on longitudinal momentum fractions from that on the external (large) momentum scale Q^, 
the latter being controlled by 



d_ ^ 



with solutions 



Xi 1 



~ <J) It- 



In 



A2 



This allows us to write the full nucleon distribution amplitude in the form 



(5) 



(6) 



$ (^X^, Q^) ~ X1X2X3 Bn (/^^) {Xi) lu 



n=0 



A2 



(7) 



where $n(xj) are appropriate but not tabulated polynomials, and the expansion coefficients 
Bn encode the nonperturbative input of the bound-states dynamics at the factorization scale 
/i. 
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From the factorized form of (xj, Q^) in Eq. (|^), it follows that the evolution equation 
for the x-dependence reduces to the characteristic equation 



X1X2X3 



3Cf _ 



$ (xi) = — / [dy] V (xi, Vi) $ (y, 



(8) 



with Cb = (iVc + 1) /2A^c. 

To proceed, it is convenient to consider the kernel V (xj, Ui) as being an operator expanded 
on the polynomial basis 1x5^X3) = \kl) (recall that because of momentum conservation, 
only two out of three Xj variables are linearly independent), i.e., to write 



and convert Eq. 



/ [dy]V{x,,y,) 
Jo 

into the algebraic equation 

'^n 'Jn (-^i 



(9) 



3 Of ,Cb V 

2 Po /3o 2w (Xi 



(10) 



where w (xj) = X1X2X3 = xi (1 — xi — X3) X3 is the weight function of the orthogonal basis. 
In this way, the action of the operator V can be completely determined by a matrix, namely: 



V\kl) 
2w (xj) 



i+j<M 



ij,kl 



(11) 



The corresponding eigenvalues are then determined by the roots ?7„ of the characteristic 
polynomial that diagonalizes the matrix U: 



V $„ (Xi) = -Tjn W (Xi) $„ (Xi) 

SO that the anomalous dimensions are given by 

Po \2 

with the orthogonalization prescription 



1 ^ ^ 1 

[dx] W (Xi) (Xi) $„ (Xi) = — 6mn , 



(12) 



(13) 



(14) 



where Nm are appropriate normalization constants (see Table |l|). 

Within the basis \kl), the matrix U can be diagonalized to provide eigenfunctions, which 
are polynomials of degree M = k + I = 0,1,2,3..., with M + 1 eigenfunctions for each M. 
This will be done within a basis of symmetrized Appell polynomials defined by 



k+l<m+n 

E 

kJ=0 



^kl 



kl) 



(15) 
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where + refers to the case m > n and — to the case m < n, and the Appell polynomials are 
defined in terms of special hypergeometric functions, according to 

^™„(xi, xs) = ^(f^Hs, 2, 2; xi, X3) • (16) 

These polyomials constitute an orthogonal set on the triangle T = T{xi, X3) with xi > 0, 
X3 > 0, Xi + X3 < 1 and provide a very suitable basis to solve the eigenvalue equation for the 
nucleon because V (i) is blockdiagonal for different polynomial orders, and (ii) commutes 
with the permutation operator P13 = [321], i.e., [P13, V] = 0. Hence V becomes, in addition, 
blockdiagonal within each sector of permutation-symmetry eigenfunctions at fixed order M. 
As a result, the kernel V can be analytically diagonalized up to order seven, providing a 
total of n^axiM) = |(M+l)(M + 2) = 36 eigenvectors and eigenvalues (alias eigenfunctions 
and anomalous dimensions). Beyond that order, the evolution equation has to be solved 
numerically because the roots of the characteristic polynomial of matrices with rank four 
cannot be determined analytically. 

From the practical point of view, such a large set of eigenfunctions surely exceeds the 
number of theoretical constraints on the nonperturbative input coefficients S„ that can be 
derived from QCD sum rules ||1| J2^J23[| , lattice simulations p4[, or by fitting experimental 



data [^. However, yet much larger orders are necessary in order to extract the asymptotic 
behavior of the anomalous- dimensions spectrum. 

For this purpose, let us express all nucleon eigenfunctions as linear combinations of 
symmetrized Appell polynomials of the same order M 

m+n=M 

^kixi)= 4n^^n(5,2,2;xi,a;3) (17) 

m,n=0 

and rearrange jF^n (cf. Eq. (|T5D), which belongs to a definite symmetry class Sn = ±1 
within order M, in the form of an (arbitrary) vector to read 

^mn {Xi,X3)\ >:Fg{Xi,X3). (18) 

Then Hilbert-Schmidt orthogonalization yields a basis 

\^',) = T.ZUkl) (19) 



with 



so that 



/ [dx] W (Xi) Tg^g, OC 5gqf , (20) 




V\J^') 1 

' o E Zl,U,^M\^J)■ (21) 



q'_ _ _ 
\ 9 

i,j,k,l 



Note that the construction of polynomials depending on two variables via the Hilbert- 
Schmidt method has no unique solution, but depends on the order in which the orthog- 
onalization is performed. Since beyond order M = 3, neither the eigenvalues nor the 
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normalization factors are rational numbers, one has to find which representation is more 
convenient for calculations. 

The last step in determining the eigenf unctions and eigenvalues of V is to define the 
matrix 



Mq'q 



[dx] W (Xi, (1 - Xi - Xa), Xg) P, (Xi, X3) V (Xi, Xg) P (Xi, X3) 



and calculate the roots of the characteristic polynomial 

V{r]) = det [Mg>g - r]Ig'g] 



(22) 



(23) 



Consequently, in terms of the eigenvectors rUg = (m^, . . . m^'^ of A^^'g, the eigenfunctions 
of the evolution equation are given by 



^g{xi,xs) oc ^m^'j^^, (xi,X3 

k,l 



(24) 



For every order M, there are M + 1 eigenfunctions of the same order with an excess of 
symmetric terms by one for even orders. A compedium of the results up to order M = 4, 
yielding a total of 15 eigenfunctions and associated anomalous dimensions, is given in Table |l|. 
The precision of orthogonality is at least 10^*^. 

It turns out that the eigenfunctions {$n} of the nucleon evolution equation satisfy a 
commutative algebra subject to the triangular condition \0{k) — 0{l)\ < 0{m) < 0{k) + 
0{l): 



(Xi) (X,) = J2 F^l^rn (Xi) 
m=0 

with structure coefficients given by 

Ffc7 = f [dx] xiX3(l - xi - X3)$m (xi) (xi) (xi) 
Jo 

0{k) being defined by 



(25) 



(26) 



0{k) 



r 

1 
2 
3 
4 



A; = 
l<k<2 
3<k<5 
Q<k<9 
k = 10,11 



(27) 



Note that the structure coefficients are symmetric, i.e., = FIJ}. Furthermore, Fj^/^ = 
The utility of this algebra derives from the fact that once the structure coefficients have 
been computed, they can be used to express any function /(xi, X3) in terms of the nucleon 
eigenfunctions. The values of F^™ up to 0{k) = 11 are tabulated in |jl8|. Unfortunately, it is 
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Table 1. Orthogonal eigenfunctions (xj) = J2kl ^kl ^1^3 of the nucleon evolution equation 

up to order M = 4 in terms of the coefficient matrix a^^ (a^^ = Sn afi^ with n fixed) and the 
corresponding anomalous dimensions 7„(M) defined in the text. The numerical results for n > 12 
have been obtained with a much higher numerical accuracy than shown in this table. 



M 



In 



Vn 






"oo 


-1 


120 




1 


2 
3 


1260 







1 


420 




-2 


6 
? 

2 


756 




2 


34020 







1944 




2 


-(-79+v^) 

24 
79+V97 

, 24 ^ 


4620 (4r>r,+ ii v 97) 
97 

4620 (485-11 V9f) 

, 97 

27720 (33607-247 ViStJT) 


1 
1 





120 
379+^1801 
, 120 


4801 

27720 (33607+247 
, 4801 


v^SoT) 


-(-256+^1081) 


196560 (7567-13 


/lOST) 





153.37061 


60 

256+</l08T 
60 

3.23876 


, 1081 
196660 (7567+13 i 

1081 
1 


/lOST) 


3.94397 


1 




332.500864 


4.81727 


1 




-137.11538 



9 

10 

11 
12 
13 
14 



1 
-1 

1 

1 
-1 

1 

1 

1 

-1 

-1 

-1 

-1 
1 
1 
1 



27 

115-v/97 

162 
115+VW 

162 

559-v/4S0l 

810 

559+V^Sor 

810 

346-\/l08T 

405 

346+v^i)8T 

405 

0.70204 
0.80651 
0.93589 



n 
~0 

1 
2 
3 
4 
5 

6 
7 
8 
9 
10 

11 
12 
13 
14 



-'lo 



-^20 



"30 


"21 






































-5(l7+v^) 


-5 (31- v^) 


16 

-5 (17-V97) 


-5 (31+V97) 
8 


16 


601+v''480T 


59-V480T 


264 


44 


eoi-v^Sn 


59+v^8f)T 


264 


44 


379+\/T()8T 


61-^1081 


48 


8 


379- vTosT 


6i+VT()8T 


48 


8 



^40 



^.31 



-7 

-6 
-6 

1 

1 

1 

1 



-1380.33552 
-2992.50778 
1234.03849 






8 

_ 4 

l| 
3 

41+^97 

4 

41-^97 
4 

-3 
-3 
-5 



5232.86956 
9240.51876 
-1843.05428 






4 

14 

3 (31-v^) 

4 

3(31+^/97) 
4 






5006.42414 
17166.06044 
-12981.41464 



-8063.85349 
-11695.76593 
-587.61051 



-9178.44426 
-31471.11081 
23799.26017 
















-(159+v'lOST) 
40 

-(159- vTosT) 
40 

4345.63139 
5068.49438 
1382.85660 
















-(61- v'108l) 
8 

-(ei+vTosT) 

4926.80699 
19489.65169 
-10302.90296 


















8503.27454 
23962.91822 
-26992.71442 



7 



J I I I I I I I I I I I I I I I I I I I I I I I I I I L 



10 20 30 

M 



Figure 1. Illustration of the spectrum of the anomalous dimensions up to order M = 30, comprising 
symmetric (full circles) and antisymmetric (open circles) contributions. 

not that easy to identify the proper group (and physical symmetry) underlying this algebra. 
Work in this direction is still in progress. 

Next, we present the results obtained for the anomalous-dimensions spectrum. An im- 
pression of the detailed structure of the spectrum is provided by fig. |l|. Both symmetry 
classes under the permutation P13 are shown: open circles stand for values belonging to 
S'„ = — 1 (antisymmetric sector) and black dots for those belonging to S'„ = 1 (symmetric 
sector). As the order M increases, the degeneracy of the spectrum also increases because 
more and more operators with the same quantum numbers contribute and the density of 
eigenvalues becomes very high. Because all 7n are positive fractional numbers increasing 
with the counting index n, higher terms in the eigenfunctions decomposition of the nucleon 
distribution amplitude are gradually suppressed (cf. Eq. |^). 

At very large order a different picture for the large-order behavior of the spectrum 
of anomalous dimensions develops, namely one of logarithmic rise. Indeed, while the 
low-order spectrum with M < 3 can be reproduced by the empirical power-law 

7„(M) = 0.37M°-^^^, the inclusion of large orders seems to be better described by a loga- 
rithm. Extending the calculation up to a maximum order of M = 400, we finally obtain the 
spectrum displayed in fig. ^ 

The anomalous- dimensions spectrum can be reproduced by the logarithmic fit 

7„(M) = c + rfln(M + 6) . (28) 

Up to order 150 both sectors of eigenvalues are included, i.e., Sn = ±1- Beyond that 
order, for reasons of technical convenience, only the antisymmetric ones have been taken 
into account. However, since the spacings become very dense asymptotically, this poses no 
restrictions on the validity of the calculation. The upper envelope of the spectrum is best 
described by the following values of the parameters with their errors: h = 1.90989 ±0.00676, 
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Figure 2. Spectrum of the anomalous dimensions of trilinear twist-3 quark operators up to order 
M = 400. The solid lines (upper and lower envelopes of the spectrum) represent logarithmic fits up 
to the maximum considered order 400, taking into consideration all orders above 10. The dashed line gives 
for comparison a previous logarithmic fit (see text) which takes into account all orders up to 150. 



c = -0.637947 ± 0.000634, and d = 0.88822 ± 0.000119. For the lower envelope, the 
corresponding values are b = 3.006±0.483, c = -0.3954±0.0290, and d = 0.59691±0.00545. 
The spacing of eigenvalues at very large order is reproduced by the values b = — 0.027±0.728, 
c = -0.2460 ± 0.0248, and d = 0.291883 ± 0.00475. It is worth remarking that the previous 
logarithmic fit 7„(M) = [logiQ(2.13M + 1.4)]^'^^, given in |12|, which takes into account all 
orders up to 150, deviates from that in Eq. (|28|) only by an amount less than 5% at order 
400. Hence, we conclude that the calculated spectrum shows already asymptotic behavior, 
rendering the inclusion of still higher orders superfluous. In particular, as one observes 
from fig. ^, the exponent d for the upper envelope shows scaling behavior for M > 200, 
approaching fast the value 0.8882. Physically, the logarithmic rise of the spectrum is due to 
the enhanced emission of soft gluons, reflecting the fact that the probability for finding bare 



quarks decreases [26 



Let us now turn to the recent work by Braun et al. |^|. These authors claim that 
the lowest two eigenvalues of the spectrum decouple, being separated from the others by 
a gap, which they interpret as resulting from the formation of a diquark state. Let us see 
whether our exact large-order calculation confirms their finding. Fig. |^ shows the gap at 
each particular order between the lowest two anomalous dimensions and the 3rd one as a 
function of the order M (marked by black dots) in comparison with the gap between the 
7th and 8th anomalous dimensions (denoted by crosses), which we found to be the largest 
gap after the 3rd one. One observes that indeed beyond order 50, the gap of the two lowest 
anomalous dimensions remains constant (up to the maximum considered order 150), whereas 
the gap of the higher ones decreases logarithmically like 0.267133/ In M (solid line), where 
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Figure 3. Determination of the exponent d for the upper branch of the anomalous-dimensions spectrum. 



the first four sliown values were fitted. Note that beyond order 200 only the antisymmetric 
values were used, so that the calculated gaps in this sector lie slightly higher than the gaps 
between the combined set of values of anomalous dimensions. This behavior can indeed be 
interpreted as supporting the results of [pl|], though our exact calculation gives for that gap 



a somewhat higher value (cf. fig. ^ than their approximate calculation. 

In summary, we have provided convincing evidence that the large-order behavior of the 
anomalous-dimensions spectrum of trilinear twist-3 quark operators with total derivatives 
tends to grow logarithmically with order. The exponent for the upper envelope of the 
spectrum was determined with high precision and found to have the asymptotic value 0.8882. 
The lowest two levels of the spectrum seem to be separated from higher ones by a gap - 
as found by Braun et al. pl| - possibly indicating the formation of a binary quark system 



(diquark cluster) inside the nucleon. Yet, scepticism remains whether this finding will survive 
the inclusion to the kernel of the nucleon evolution equation of higher-order contributions 
which comprise gluon self-interactions. We have presented explicit results for a set of 15 
orthonormalized eigenfunctions of the nucleon evolution equation, employing a symmetrized 
basis of Appell polynomials, whereas still higher states can be readily constructed by the 
means developed and used in this work. These results suggest that the objections raised in 
|2T| against the use of Appell polynomials in solving the nucleon evolution equation are not 
really justified. 
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Figure 4. Gaps of eigenvalues in the spectrum of anomalous dimensions. 
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